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Abstract 

Using Leaver's continue fraction and time domain method, we investigate the wave dynamics 
of phantom scalar perturbation in the background of Schwarzschild black hole. We find that the 
presence of the negative kinetic energy terms modifies the standard results in quasinormal spectrums 
and late-time behaviors of the scalar perturbations. The phantom scalar perturbation in the late- 
time evolution will grow with an exponential rate. 
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Many observations confirm that our universe is undergoing an accelerated expansion. In order to explain this 
observed phenomena, the universe is supposed to be filled with dark energy within the framework of Einstein's 
general relativity. Dark energy is an exotic energy component with negative pressure and constitutes about 
72% of present total cosmic energy. The simplest explanation for dark energy is the cosmological constant [l| , 
which is a term that can be added to Einstein's equations. This term acts like a perfect fluid with an equation 
of state ujx = —1, and the energy density is associated with quantum vacuum. Although this interpretation 
is consistent with observational data, it suffers the coincidence problem, namely, "why are the vacuum and 
matter energy densities of precisely the same order today?". Therefore the dynamical scalar fields, such as 
quintessence 3] i k-essence [3| and phantom field J] , have been put forth as an alternative of dark energy. 

The phantom field model is an interesting candidate for dark energy since it has some peculiar properties. 
This field has a negative kinetic energy and so that the null energy condition is violated and the equation 
of state LOx of dark energy less than —1. The super negative equation of state is favored by recent precise 
observational data involving CMB, Hubble Space Telescope, type la Supernova, and 2dF data sets 51. The 
dynamical properties of the phantom field in the cosmology has been investigated in the last years a III ISl, IS , 
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13[ . It shows that the energy density increases with the time and approaches to infinity in a finite 
time 6]. In other words, the universe dominated by phantom energy will blow up incessantly and arrive at 
a big rip finally, which is a future singularity with a strong exclusive force so that anything in the universe 
including the large galaxies will be torn up. Recently, many efforts have been working to avoid the big rip 



14| . It has argued that if one considers the effects from loop quantum gravity, this future singularity will be 



disappeared in the universe 
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It is of interest to extend the study of dynamical evolution of the phantom field to black hole spacetime 
in various gravity theories, since this could help us to obtain the connection between dark energy and black 
hole. The dynamical evolution of usual scalar field perturbation has been studied for the last few decades 
(for a review, see [l^ Q, lill). It is well known that a static observer outside a black hole can indicate three 
successive stages of the wave evolution. The first one is that the exact shape of the wave front depends on 
the initial pulse. This stage is followed by a quasinormal ringing, which describes the damped oscillations 
under perturbations in the surrounding geometry of a black hole with frequencies and damping times of the 
oscillations entirely fixed by the black hole parameters. It is widely believed that the quasinormal modes 
carry characteristic fingerprint of a black hole and can offer a direct way to identify the black hole existence. 



Detection of these quasinormal modes is expected to be realized through gravitational wave observation in the 
near future 
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2C|. Apart from the potential astrophysical interest, quasinormal modes could also serve as a 



tool to test ground of fundamental physics. It has been argued that the study of QNM can help us get deeper 
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dS/CFT 



23| correspondences, loop quantum gravity 



241 and also 



understandings of the AdS/CFT 
the phase transition of black holes 25|. In this letter, we treat the phantom field as an external perturbation 
and examine whether there exists some new feature in the dynamical evolution of phantom scalar field in a 
black hole spacetime. 

In the curve spacetime, the action of the phantom scalar field with the negative kinetic energy term is 

R 1 



S 



WttG 



(1) 



Here we take metric signature (H ). The usual "Mexican hat" symmetry breaking potential has the form 

1 



(2) 



where fj, is the mass of the scalar field and k is the coupling constant. In general, the presence of the phantom 
scalar field will change the structure of the black hole spacetime 26|. Here we just treat it as an external 



perturbation and suppose it does not affect the metric of the background. Meanwhile we only consider the 
case K = for conveniently, namely, the potential has the form V{tp) = —^fi^t/j'^. The wave dynamics of usual 
scalar perturbations with this type of the potential form has been extensively investigated in the various black 
holes spacetime 
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28|. 



Varying the action with respect to ip, we obtain the wave equation for phantom scalar field in the curve 
spacetime 
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(3) 



Comparing the Klein-Gordon equation of usual massive scalar field, we find the unique difference in equation 
([3]) is that the sign of the mass term /i^ is negative, which will yield the peculiar dynamical evolution of the 
phantom perturbation in the black hole spacetime. 

Let us now to consider the case of a Schwarzschild black hole spacetime, whose metric in the standard 
coordinate can be described by 



ds^ ^ (1 _ ^)df - (1 - —y'dr' - r^de' + sin' Odcf"). 



(4) 



Separating ip — e ^'^'^R{r)Yi„i(9,(p)/r, we can obtain the radial equation for the scalar perturbation in the 



Schwarzschild black hole spacetime 



d^Rjr) 
drl 



+ [uj"^ - V{r)]R{r) = 0, 



(5) 



where r, is the tortoise coordinate (which is defined by dr 



r-2M '^^) ^^"^ ^^'^ effective potential V{r) reads 

, , / 2M\ f 1(1 + 1) 2M .\ 
nr)^(l- — ][^^ + — -f^')- (6) 



Obviously, as /i = the radial equation of phantom scalar field reduce to that of usual massless scalar field. 
In the case fj, 0, the effective potential V(r) vanishes at the event horizon and approaches to a negative 
constant — /i^ at the spatial infinity. This is different from that of usual massive scalar perturbation. Moreover, 
from figure (1), one can see that as the mass /i increases the effective potential V{r) for the phantom scalar 
perturbation decreases and for the usual one increases. This implies the wave dynamics of the phantom scalar 
perturbation possesses some different properties from that of the usual scalar perturbation. 





FIG. 1: Variety of the effective potential V{r) with r for different jj,. The left for phantom scalar field and the right 
for normal scalar field. Here I = and M — 0.5. 



We are now in a position to apply the continue fraction method fSCrl and calculate the fundamental quasi- 
normal modes (n = 0) of phantom scalar perturbation in the Schwarzschild black hole. From equation ([5]), 
we know that the boundary conditions on the wave function R{r) have the form 

(r-I)-*", r^I, 



(7) 



r 2x e 



oo, 



where we set 2M — 1 and x = ^/t^'^ + /i^. A solution to Eq.(l5]) that has the desired behavior at the boundary 
can be written as 



R{r) 



71=0 



r~ 1 



(8) 



The sequence of the expansion coefficients {a„ : n = 1, 2, ....} is determined by a three-term recurrence 
relation staring with qq — 1: 



aoQi + /3oao = 0, 

Q!„a„+i + /3„a„ + 7„a„_i =0, n = I, 2, 

The recurrence coefficient a„, /?„ and 7„ are given in terms of n and the black hole parameters by 

an = n^ + {Co + l)n + Co, 

(}n^-2n^ + {Ci+2)n + C3, 

7„ = + (Ca - 3)n + C4 - C2 + 2, 



(9) 



(10) 



and the intermediate constants C„ are defined by 



Co = 1 - 2iw, 



Ci = -4 + i(4tj + 3x) + 



C2 =3-i(2L^ + x) 



^3 = 
C4 = 



c^ + X 



X 
2x 



, ,0 i(u) + 3x) 



-Z(/ + f)-f, 



+ i 



(11) 



If the boundary condition ([7]) is satisfied and the series in ([5]) converge for the given I, the frequency w is a 
root of the continued fraction equation 



/3« 



a„-i7ri 



. a 

This equation is impossible to be solved analytically. But we can rely on the numerical calculation to obtain 
the quasinormal frequencies for phantom scalar perturbations in the Schwarzschild black hole spacetime. In 



Qti + 177i + 2 



(12) 
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I = 1 


I = 2 





0.220910-0.2097911 


0.585872-0.195321 


0.967288-0.1935181 


0.1 


0.223315-0.2275011 


0.579967-0.2006251 


0.962955-0.1955991 


0.2 


0.228479-0.2693821 


0.562925-0.2170491 


0.949962-0.2019951 


0.3 


0.232013-0.3212981 


0.537016-0.2456961 


0.928354-0.2131751 



TABLE I: The fundamental (n = 0) quasinormal frequencies of phantom scalar field In the Schwarzschild black hole 
spacetime. 



table I, we list the fundamental quasinormal frequencies of phantom scalar perturbation field for fixed I = 0, 
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FIG. 2: Change of the real part quasinormal frequencies (n = 0) for the phantom scalar perturbations with the 
parameter /i for fixed 
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FIG. 3: Change of the imaginary part quasinormal frequencies (n = 0) for the phantom scalar perturbations with the 
parameter jj, for fixed I. 
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I — 1 and ^ = 2 in the Schwarzschild black hole spacetime. From the table I and figures (1) and (2), we find 
that with the increase of the mass /i the real parts increase for / ~ and decrease for Z = 1 and 1 = 2. The 
absolute value of imaginary parts for all I increase, which can be explained by that the bigger fj, leads to lower 
peak of the potential and thus it is easier for the wave to be absorbed into the black hole. The dependence of 
quasinormal modes on the mass fi is different from that of the usual scalar field since their effective potentials 
are quite a different in the black hole spacetime. 

Now we will extend to study the late-time behavior of phantom scalar perturbations in the Schwarzschild 
black hole spacetime by using of the time domain method [sij]. 



Adopting the null coordinates u = t — and v = t + r^, the wave equation 



(13) 



can be recast as 



dudv 



+ T/(r)V' = 0, 



(14) 



The two-dimensional wave equation (|14p can be integrated numerically, using for example the finite difference 
method suggested in 



31[. Using Taylor's theorem, it is discretized as 



Wn ^i^E + Vw - Ws- dudvV{ ) h O(e^) = 0, 



(15) 



where the points iV, S, E and W from a null rectangle with relative positions as: N: {u + Su^v + 5v), 
W: {u + 5u,v), E: {u,v + Sv) and S: {u,v). The parameter e is an overall grid scalar factor, so that 
Su Sv ^ e. Considering that the behavior of the wave function is not sensitive to the choice of initial data, 
we set ip{u, V ~ vq) = and use a Gaussian pulse as an initial perturbation, centered on Vc and with width a 
on u = as 



'ip{u — uq, v) = exp[— - 



2a2 



(16) 



Here we confirm numerically the peculiar properties in late-time evolution of phantom scalar perturbations in 
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FIG. 4: The late-time behaviors of the phantom scalar perturbations for fixed I — 0, the left and the right are for 
fi = 0.01 and 0.02, respectively. The dashed line denotes the function ln?/j ~ afit — 4/ — /3. We set the numerical 
constants a = 0.27, /3 — 3.93 in the left figure and a = 0.36, /3 = 3.92 in the right one. The constants in the Gauss 
pulse (|16p Vc — 10 and (7 = 3. 



the Schwarzschild black hole spacetime. In figures (4), (5) and (6), we plot the evolutions of the phantom scalar 
perturbations for fixed I = 0, I = I and I = 2, respectively. Unlike the usual massive scalar perturbations, 
the phantom scalar field after undergoing the quasinormal modes stage does not decay but grow with the 
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FIG. 5: The late-time behaviors of the phantom scalar perturbations for fixed / — 1, the left and the right are for 
fi = 0.01 and 0.02, respectively. The dashed line denotes the function In?/; ~ afit — 4/ — /3. We set the numerical 
constants a = 0.20, f3 — 4.32 in the left figure and a = 0.32, /3 = 3.84 in the right one. The constants in the Gauss 
pulse (|16p Vc — 10 and (7 = 3. 
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FIG. 6: The late-time behaviors of the phantom scalar perturbations for fixed 1 — 2, the left and the right are for 
H = 0.01 and 0.02, respectively. The dashed line denotes the function In-i/i ~ afit — 4/ — /3. We set the numerical 
constants q = 0.18, fH — 5.55 in the left figure and a — 0.28, /3 = 4.30 in the right one. The constants in the Gauss 
pulse (|16|l Vc — 10 and (7 = 3. 



time. The asymptotic behaviors of wave function for the phantom scalar field in the Schwarzschild black hole 
spacetime can be fitted best by 



■0 ~ e' 



ant-Al-0 



(17) 



where a and /? are two numerical constant. It means that the phantom scalar perturbation grows with 
exponential rate in the late-time evolution. This behavior can be attributed to that the effective potential is 
negative at the spatial infinity shown in fig.(l), which implies that the wave outside the black hole gains energy 
from the spacetime. For the larger /i, the scalar perturbation grows more faster since the larger /i corresponds 
to more negative potential. It is quite different from that of the usual massive scalar perturbations in a black 
hole spacetime which decay with the oscillatory inverse power-law behavior t^^ sin fit Moreover, the 

growing modes of the late-time tails caused by the negative effective potential was also observed in [29 1 . 



The exponential growth of the phantom scalar perturbation in the late-time evolution also means that it 
is unstable in the black hole spacetime and its energy density will increase with the time, which is similar to 
that of in the Einstein cosmology. Our result also agrees with that obtained in the phantom energy accretion 
of black hole [32| , where the black hole mass will be decreased and the energy density of the phantom will 



increase. 

In summary we examine the wave dynamics of the phantom scalar perturbation in the background of 
Schwarzschild black hole. Our results show that due to the presence of the negative kinetic energy, the 
properties of the wave dynamics of phantom scalar perturbation are different from that of the usual massive 
scalar field. In the late-time evolution the phantom field does not decay but instead grows with an exponential 
rate. It would be of interest to generalize our study to other black hole spacetimes, such as rotating black 
hole and stringy black holes etc. Work in this direction will be reported in the future. 
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